Evidence for strong Coulomb correlations in metallic phase of vanadium dioxide 
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The influence of Coulomb correlation on magnetic and spectral properties in metallic rutile phase 
of vanadium dioxide is studied by state of the art LDA+DMFT method. Calculation results in 
strongly correlated metallic state with an effective mass renormalization m*/m ~ 2. Uniform 
magnetic susceptibility shows Curie- Weiss temperature dependence with effective magnetic moment, 
P^tifT^ ~ l-54/is, in a good agreement with the experimental value p^jj = 1.53/is that is close to 
ideal value for V'*"'' ion with the spin S — 1/2, Peff = 1.73/iB. Calculated spectral function shows 
well developed Hubbard bands observabale in the recent experimental photoemission spectra. We 
conclude that VO2 in rutile phase is strongly correlated metal with local magnetic moments formed 
by vanadium d-electrons. 



PACS numbers: 

Transition metal compounds and particularly vana- 
dium oxides attract a lot of interests because of rich va- 
riety of structural modifications and properties^. Vana- 
dium dioxide is not an exception and is well known due 
to metal-insulator transition at 340 K— . The value of the 
energy gap (Ag=0.6 eV) in insulating phase and vicin- 
ity of the transition and room temperatures allows one 
to use this material as "wisdom" windows and in an- 
other termochromic applications^. This compound has 
been extensively investigated experimentally and theo- 
retically over the half of century. The phase diagrani^ 
in the temperature-pressure (or chemical doping) coordi- 
nates consists of the paramagnetic metal at high temper- 
ature and three insulating states below Tc- Below transi- 
tion temperature there are three insulating phases, called 
Ml, T and M2 at low, intermediate and large values of 
doping, respectively. The Mi and M2 phases crystallize 
to monoclinic structure which can be obtained from high- 
temperature rutile by distortion and doubling of the cell. 
In the former structure equivalent vanadium atoms are 
paired and off-axis shifted forming zigzags while in the 
later one there are two kinds of vanadium chains, one 
of which is paired and another one forms zigzags with- 
out pairing. In both phases the conductivity drops by 
factor of several orders in magnitude comparing with ru- 
tile phase with the estimation of the energy gap about 
0.6 eV— . Above the transition temperature VO2 crystal- 
lizes to rutile structure at all values of doping and can 
be regarded as two vanadium chains along cr axis (see 
Fig.P). 

Most of earlier theoretical studies for correlation effects 
in vanadium dioxide were concentrated on the insulating 
low-temperature phases where band structure calcula- 
tions have difficulties to reproduce correct insulating type 
of electronic structure^"—. Calculation for rutile crystal 
structure resulted in metallic bands in agreement with 
experiment, and hence, it was assumed that taking into 
account Coulomb correlations is not needed here. How- 
ever magnetic and spectral measurements for VO2 show 
that correlation strength could be essential for metallic 
rutile phased. The magnetic experiments for metallic 



rutile phase show a quite large temperature dependent 
magnetic susceptibility that indicates formation of local 
magnetic moment in the partially filled d-band. Recent 
photoemission datai^ display the quasi-particle peak at 
the Fermi level and high energy structure associated with 
the lower Hubbard band about -1 eV. The effective elec- 
tron mass, m*/m, extracted from optics^ is about 4.3 
which is larger then the value observed for Sr^jCai-ajVOa 
characterized as strongly correlated paramagnetic met- 
alsiS. Recent LDA-I-DMFT investigationsii of vanadium 
oxide have paid an attention to the spectral properties of 
the compound and have showed a good agreement with 
photoemission and optical data. 

In the present work we have used state of the art 
LDA+DMFT methodic to study the physics of VO2 in 
metallic rutile phase with the special focus on the mag- 
netic properties. Temperature dependence of the cal- 
culated magnetic susceptibility xC^) reproduces with a 
good accuracy Curie- Weiss law xcw{T) = plff/3{T—Q) 
and is in an agreement with experimental data. Calcu- 
lations results for the spectral function are typical for 
strongly correlated metal with well developed Hubbard 
bands in addition to quasi-particle metallic peak at Fermi 
level and they agree well with earlier studiesii. Calcu- 
lated value of effective electron mass enhancement factor 
m* /m « 2 proves a relatively large strength of correla- 
tion effects for metallic phase of VO2. Therefore, vana- 
dium dioxide in rutile phase is the strongly correlated 
metal with local magnetic moments. 

At high temperatures, vanadium dioxide crystallizes as 
rutile shown on the Fig. [TJ It has the PA2/mnm space 
group and details of lattice constants and atomic posi- 
tions can be found in Ref. [Tsl . Each vanadium atom is 
surrounded by a oxygen's octahedron resulting in split- 
ting of d level to triply degenerate 2g and doubly degen- 
erate states. Additional tetragonal distortion present 
in the structure leads to further lifting of degeneracy of 
2g level to aig and tt states (d|| and (i,r according to 
Gooenough's notations^). In the local coordinate basis 
shown on the Fig. [1] by (double) primed axes for vana- 
dium atoms, aig state corresponds to the d^'y' orbital. 



FIG. 1: (Color online) Crystal structure of rutile phase. Vana- 
dium and oxygen atoms are denoted by large (red) and small 
(cyan) balls. {aR,bR, cr) axe rutile crystallographic axes and 
coincide with the Cartesian one. Local coordinate axes on the 
vanadium atoms are shown by (double) primed set. 



TT states correspond to d^'z' and dyiz' and to dz>2 and 
dx'2-yi2 orbitals. The local z' axis is chosen to be pointed 
to apical oxygen atoms and local x' and y' are chosen to 
be pointed to the planar oxygens. This local coordinate 
system can be obtained by the rotation to the Euler's 
angles (7r/4, — 7r/2, 7r/4). The double primed local coor- 
dinate axes are obtained by a proper symmetry operation 
for corresponding vanadium atoms. 

We used TB-LMTO methodic to calculate electronic 
structure of VO2. The results obtained are presented on 
the Fig. [5] and have a good agreement with earlier stud- 
iesi. Total and partial (per atom) density of states are 
shown on the top panel (for details of color-coding see 
caption). One can clearly see that the states crossing 
Fermi level arc of vanadium 2g symmetry mainly and 
spread from -0.5 eV to 2 eV. They are well separated by 
a gap of 1.4 eV from occupied bonding states of mixed 
O 2p and V character. The anti-bonding combination 
is located above 2 eV and is also separated by a tiny 
gap from the 2g bands. The orbitally (symmetry) de- 
composed partial DOScs for vanadium atom in the local 
coordinate system mentioned above are presented on the 
bottom panel of Fig. [51 aig state is more narrow than tt 
one and has a one-dimensional-like shape. 

We constructed low-energy Hamiltonian for V 2g states 
using the A^MTO method-^. It has a dimension of 
6 and will be used in LDA+DMFT calculation. The 
eigenvalues of the iVMTO Hamiltonian coincide with the 
LMTO's one by construction in the energy interval close 
to the Fermi level and shape of the density of states re- 
mains the same (not presented). Hopping integrals for 
the Wannier function of aig, tti and 1:2 idx'y',dx'z' and 
dy'z' in local notations) symmetries are presented in the 
Table |I] for the vanadium atoms in chain (along cr axis) 
and between chains (upper and lower parts of table, re- 
spectively). One can immediately see that aig electron 




FIG. 2: (Color online) LDA density of states for the rutile 
phase of VO2. Top panel shows total, V 3d and O 2p DOSes 
by (black) solid, (green) dashed and (red) dot-dashed lines, 
respectively. Bottom panel shows partial DOSes of different 
symmetry for vanadium atom (see legends for notations) 



can easily jump to the same orbital on the nearest neigh- 
bor atom and it is the largest (intra chain) hopping. The 
probability to hop to another chains is much smaller for 
aig electron. Overall structure of aig hoppings assumes 
one-dimensional-like physics where the largest transfer 
integral defines energy splitting between peaks and the 
rest is responsible for asymmetry of density of states and 
its broadening (see aig partial DOS on the bottom panel 
of Fig. [2]). The structure of the tt hoppings is more 
uniform in space that produces more three dimensional 
shape of density of states. 

To take into account electron-electron correlation 
in the partially filled d-band we use LDA+DMFT 
methodic. It combines the material specific aspect in 
LDA and accurate treatment of local Coulomb intrac- 
tion in BMFT-^^. The later is achieved via mapping 
of the lattice problem with many degrees of freedom 
to the quantum impurity embedded to the time depen- 
dent self-consistent bath (for details about LDA-I-DMFT 
method and implementations see e.g. Refs fisl) . For the 
solution of the quantum impurity problem the modern 
multi-orbital Hirsch-Fye quantum Monte-Carlo method 
has been used^^. In our LDA+DMFT calculations we use 
the value of screened Coulomb interaction, U =4 eV, and 
the value of Hund's exchange, J=0.68 cV. In the quan- 
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TABLE L Hopping integrals for the rutile phase (in meV). 
The distance between vanadium atoms is in second column. 
The largest hoppings to another distances are smaller than 
25 meV and not shown here. 



turn Monte Carlo calculations we keep Ar = /3/L=0.25 
for all temperatures and of order 10*' ~ 10^ QMC steps 
to satisfy ergodicity. 

We start the presentation of the LDA+DMFT results 
from magnetic properties. Following the linear response 
idea, the uniform magnetic susceptibility can be calcu- 
lated directly by adding to the Hamiltonian an external 
magnetic field, ^bHz, and measuring magnetization of 
compound. Then, the uniform magnetic susceptibility is 



(1) 



where m{T) — 



Tl ^ 

m rn 



is a magnetization at 



given temperature. The few magnetic fields of order 
^ 0.02-^0.1 eV were used to check and satisfy the condi- 
tion of linearity of magnetization. The inverse of calcu- 
lated uniform magnetic susceptibility and experimental 
data extracted from Ref . [l^ are shown on the Fig. [3] One 
can clearly see that these inverse quantities agree accu- 
rately to the shift along ordinate axis. This implies that 
the calculated and experimental effective magnetic mo- 
ments are in good agreement. The theoretical estimate 
obtained by the fit to the Curie- Weiss law 



Xcw{T) 



Peff 



3(T-e) 



(2) 



gives the value of the effective magnetic moment, P*//" — 
1.54/xb, in a good agreement with the experimental value 
p^Yf ~ 1.53/iB that is close to ideal value for V^"*" ion with 
the spin S = 1/2, pe// = 1.73^b. To make the fit shown 
by lines we have used the data points above 500 K where 
the linearity of is more pronounced for theoretical 
and experimental data as well. The use of complete data 
set leads to the slightly higher values of magnetic mo- 
ments^i. The theoretical Curie- Weiss temperature Q is 
almost twice larger then its experimental counterpart (- 
1600 K and -700 K correspondingly), that is connected 
to the mean-field nature of the calculations and, hence, 
absence of the collective excitations. 

We have also calculated a local spin-spin correlation 
function, {Sz{t)Sz{0)) at different temperatures shown 



Temperature (K) 

FIG. 3: (Color online) Inverse of uniform magnetic suscep- 
tibility. LDA-f DMFT results are shown by (black) squares, 
the experimental data by Zylbersztejn et ali^ are (red) cir- 
cles. Lines are result of linear fit of data above 500 K. 



on the Fig. Albeit this is a pure local quantity and 
cannot be directly compared to the experimental data, it 
gives a sound source of information about magnetic prop- 
erties of the system. At low temperatures, in the Fermi 
liquid regime, the local spin-spin correlation function be- 
haves as \T / sin{'KTT)Y' at intermediate times, r, while it 
saturates to the constant value for the insulator or metal 
with a very well developed local magnetic momen t 
At middle point, the correlator is (S'^(/3/2)S'^(0)) « T", 
with exponent a=2 for Fermi liquid and zero for the insu- 
lator. On the upper panel of the Fig. |4j from the time de- 
pendency one can see that the correlation function is not 
saturated. At the same time, the {Sz{l3/2)Sz{Q)) quan- 
tity plotted in the inset of upper panel shows the linear 
behavior with the exponent a=l different from both in- 
sulating and metallic regimes. This observation suggests 
that the formation of the local magnetic moment in the 
metallic phase is not completely finished^^. The lower 
panel of the Fig. |4] shows the analytical continuation of 
the local spin-spin correlation function to the real energy 
axis^^. The width of the peak at zero frequency is inverse 
proportional to a " lifetime " of the local magnetic moment 
and it is much larger than that of the a-iron in para- 
magnetic phased which is closer to the localized limit. 
The inset in the lower panel shows the inverse of the lo- 
cal magnetic susceptibility, xioc{T) = dT{Sz{T)Sz{0)) 
versus temperature. This quantity has a noticeable lin- 
ear behavior and by fitting it to the Curie law one ob- 
tains the value of the effective local magnetic moment, 
Pioc = 1.54^B and critical temperature. The comparison 
of the local and uniform susceptibilities gives a rough 
estimation of the degree of spacial correlations. The 
close values of the extracted effective magnetic moments, 
ptheor g^^^ pioc, assumcs that the vanadium dioxide at 
high temperatures has a localized type of magnetism. 

Now, we would like to turn our discussion to single par- 
ticle properties. The imaginary part of the self-energy for 
aig and tt orbitals are shown on the Fig. \5\ The behav- 
ior of the imaginary part for both orbitals is of Fermi- 
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FIG. 4: (Color online) The local spin-spin correlation func- 
tion for different temperatures in imaginary time domain (up- 
per panel) and its analytical continuation to the real energies 
(lower panel). Upper inset shows the dependence of this cor- 
relator at /3/2 with temperature. Lower inset is an inverse of 
the local susceptibility (for details see text). 

liquid type. At small real frequencies, in Fermi-liquid 
regime, the imaginary part is 32 (w) = —TT^ — BuP' 
that transforms to the linear behavior on the Matsubara 
axis, 3S](za;„) = —TT^ — -Bw„, observable on the figure. 
Additionally, the mass enhancement factor 

m /m =1 (3) 

as a function of temperature shows a saturation (upper 
inset of Fig. [5]) indicating the Fermi-liquid behavior— of 
moderately correlated metal with m* /m k, 1.8 and 1.7 
for aig and tt orbitals, respectively. The lower inset shows 
the extracted to zero imaginary part of self-energy, 3I](0) 
versus temperature. One can see that it goes quadrati- 
cally to zero at low temperatures once again confirming 
Fermi-liquid nature of compound. From these data one 
can conclude that the coherence temperature is about 
700 K for both orbitals. 

The comparison of the photoemission Aa.ic& (PES) 
with the calculated LDA and LDA-I-DMFT results are 
presented on the upper panel of the Fig. [Sj One can 
clearly observe a good agreement of the LDA-I-DMFT 
(color lines) and PES data (dots). There is a quasi- 
particle peak at the Fermi level and lower Hubbard band 
at higher energies. The later is an additional evidence of 



FIG. 5: (Color online) Imaginary part of self-energy for a\g 
(top) and eg (bottom) states in metallic phase. Solid (blue), 
dashed (green) and dot-dashed (red) lines present /3=30, 15 
and 10 eV"^ respectively (~390, 770 and 1160 K). Upper inset 
shows the effective mass, m*/m, for aig and eg orbitals by 
(black) circles and (dark green) squares, respectively. Lower 
insets shows the imaginary part of self-energy extracted to 
zero energy, 3E(0). Color-coding is the same as in previous 
inset. 



the correlated nature of the metallic rutile phase. Tem- 
perature dependence of the LDA+DMFT spectra shown 
by different colors is also in good qualitative agreement 
with the experiment (see Fig. 3 of the Ref. |^: the 
spectral weight is transferred from quasi-particle peak 
to higher energies and becomes smaller with decrease of 
temperature. At the same time, the LDA results shown 
by the dashed line describe the quasi-particle feature at 
the Fermi level with a large overestimation of the peak 
weight and complete absence of the high energy inco- 
herent hump. LDA-I-DMFT spectral function for lowest 
calculated temperature, T k. 390 K is presented on the 
lower panel of the Fig. [6l One can clearly see the three 
peaks structure which is typical for correlated materi- 
als. The quasi-particle peak lies at/above the Fermi level 
and low and upper Hubbard bands are located at higher 
frequencies, at -1.5 eV and ~2.3 eV. The width of the 
quasi-particle peak is reduced to ^1.1 eV comparing with 
2.1 eV value for non- interacting band width and is almost 
independent with temperature which is compatible with 
the weak temperature dependence of the quasi-particle 
weight, Z = m/m* , for both orbitals (see upper inset of 
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Fig. [5]). One should note here that the spectral functions 
for higher temperatures are similar to described above. 
The difference is that the shape of the quasi-particle peak 
becomes more complicated due to the decrease (in abso- 
lute value) of the imaginary part of self-energy at zero 
energy (see lower inset of the Fig. [5|), and hence, dis- 
playing the features of original LDA spectral function 
smoothed by high temperature at f3~10 cV~^. 



1 1 1 1 1 

1 1 


1 

* ■ 


■ ■ ■ Exp. data / i 

LDA ! A 

P = JO / A 

P = ^5 / \ 

— p = io ./ p 


\ 
\ 
\ 

\ 







-2 -1 

Energy (eV) 



2 - 



1 I 


1 1 1 1 1 
n Total 




\ CL, 

J \ ^ 










/ / 













-2 2 4 

Energy (eV) 

FIG. 6: (Color online) Top panel shows comparison of the 
photoemission data extracted from the Ref. @ and calculated 
LDA and LDA-fDMFT results for various temperatures. Bot- 
tom panel is total and orbitally resolved LDA-I-DMFT spec- 
tral functions for /3=30 eV~^ (T^390 K). 



In conclusion we have studied the influence of Coulomb 
correlation effects on magnetic and spectral properties of 
VO2 in metalhc rutile phase by LDA-I-DMFT method. 
The calculation results show typical strongly correlated 
metal behavior for self-energy with a sizable spectral 
weight transfer to Hubbard bands and electronic states 
renormalization near the Fermi level with effective mass 
m*/m w 2. The uniform and local magnetic suscepti- 
bility calculation obey Curie- Weiss temperature depen- 
dence with effective magnetic moment value close to ideal 
ionic value corresponding to configuration. A good 
agreement of calculated and measured spectral and mag- 
netic properties allows to conclude that VO2 in rutile 
phase is in strongly correlated metal state with local mag- 
netic moments formed by d-electrons. 
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